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Abstract 

s'  It  has  been  found  that  fixed  error  dead-zones 
s  defined  in  the  existing  literature  result  in 
erious  degradation  of  performance,  due  to  the  con- 
ervativeness  which  characterizes  the  determination 
f  chair  width.  In  Sfcs  piesmr  paper,  variable 
width  dead- zones  are  derived  for  the  adaptive  con¬ 
trol  of  plants  with  unmodeled  dynamics.  The  der¬ 
ivation  makes  use  of  information  available  about 
the  unmodeled  dynasties  both  a  priori  as  well  as 
during  the  adaptation  process,  so  as  to  stabilize 
the  adaptive  loop  and  at  the  same  time  overcome  the 
conservativeness  and  performance  limitations  of 
fixed-dead  zone  adaptive  or  fixed  gain  controllers,. 


1.  INTRODUCTION 


sho*\i 


Research  in  recent  years  has  shotfe  that  adap¬ 
tive  control  algorithms  which,  under  ideal  assump¬ 
tions,  have  been  proven  globally  asymptotically 
stable,  indeed  exhibit  unstable  behavior  in  cir¬ 
cumstances  under  which  chose  assumptions  are  even 
slightly  violated.  Of  the  two  instability  mech¬ 
anisms  identified  for  these  algorithms,  -commonly 
referred  to  as  "gain"  and  " phase *  instability  mech¬ 
anisms  (11#-  che  former  is  more  unavoidable  and  is 
triggered  by  the  controller  parameter 
drift  which  occurs  as  a  result  of  nonzero  output 
errors.  These  are  a  consequence  of  the  face  that, 
in  the  presence  of  uamodeled  dynamics  and/or 
(persistent)  disturbances  there  can  be  no  perfect 
(transfer  function)  matching  between  the  compensated 
plant  and  the  reference  model  over  all  frequencies, 
even  if  "sufficiency  of  excitation"  for  the  "nom¬ 
inal"  model  order  is  guaranteed. 

Perfect  matching,  on  the  other  hand,  trans laces 
into  zero  output  (tracking)  error,  under  ideal  as¬ 
sumptions,  and  has  been  the  basis  for  the  parameter 
adjustment  laws;  only  when  the  output  error  is  zero 
does  adaptation  stop.  Clearly,  then,  by  design, 
any  nonzero  output  error  is  instantaneously  attrib¬ 
uted  to  parameter  errors.  Furthermore,  there  is 
nothing  in  the  mathematics  of  the  adjustment  mech¬ 
anisms,  as  they  currently  stand,  to  pravent  gain 
drift  due  to  error  sources  other  than  parameters, 
as  for  axaaple  happens  even  in  cases  of  "exact 
modeling",  with  "sufficiency  of  excitation",  where 
convergence  to  the  "desired"  parameter  values  has 
been  achieved  momentarily;  extraneous 
disturbances  entering  at  that  point  can  cause  the 


parameters  to  drift  from  their  "desired"  values. 

Consequently  then,  the  zero  tracking  error 
requirement  must  be  suitably  relaxed.  The  ration¬ 
ale  is  that  no  parameter  adjustment  should  take 
place  when  the  output  error (s)  are  due  to  disturb¬ 
ances  and/or  unmodeled  dynamics.  This  can  be 
achieved  on  an  existing  algorithm  by  a  dead- zone 
nonlinearity,  in  the  parameter  adjustment  law, 
whose  width  depends  on  the  contribution  of  the 
disturbances  and/or  unmodeled  dynamics  to  the  out¬ 
put  error. 

The  idea  of  a  dead-zone  nonlinearity  in  the 
parameter  update  law  to  avoid  the  effect  of  disson¬ 
ances  on  adaptation  was  first  introduced  for  in¬ 
direct  adaptive  algorithms  by  Egardt  in* 1980  [2) 
and  was  later  amplified  by  Samson  [3).  Also,  in 
1962  Peterson  and  Narendra  used  a  dead-zone  non¬ 
linearity  to  prove  stability  for  a  class  of  direct 
algorithms  in  the  presence  of  bounded  disturbances 
with  no  unmodeled  dynamics  [4].  However,  the  width 
of  the  dead  zone  was  chosen  to  be  constant  and  had 
to  be  based  on  a  very  conservative  bound  so  that  it 
yialded  only  marginally  stable  systems  with  ex¬ 
tremely  poor  model  tracking  as  the  examples  in  [4] 
seem  to  suggest. 

Consequently,  obtaining  non- fixed  accurate 
bounds  for  the  disturbance  and  high-frecuency 
dynamics  contributions  to  the  output  error  is  cru¬ 
cial  to  overcoming  the  conservativeness  of  the 
dead -zone  width  which  they  define.  This  depends  on 
the  -  ability  to  translate  f requency  domain 
magnitude  bounds,  most  naturally  expressed  by 

l  -norms  into  time-domain  magnitude  bounds 
of  instantaneously  measured  quantities,  most  nat¬ 
urally  expressed  by  1*-  or,  for  our  purpbses.^L^- 
norms  which  are  much  less  conservative  than  L  - 
norms. 

This  paper  discusses  the  use  of  a  deadzone, 
whose  width  is  adjustable  on  line,  to  adaptively 
control  a  plant  with  unmodeled  dynamics,  with  the 
objective  of  maintaining  its  stability  and  minimi¬ 
zing  tne  adverse  effects  of  a  conservative  dead- 
zone  width  to  its  performance.  Due  to  space  consid¬ 
erations  we  do  not  treat  tha  case  of  output  dist*bs>- 
ances  here;  also,  the  topic  of  disturbances  addi¬ 
tionally  includes  a  fixed  disturbances  rejection 
mechanism  that  introduces  a  modification  in  the  bas¬ 
ic  structure  of  the  MRAC  system  so  is  to  merit 
separate  attention. 

Section  2  of  this  paper  contains  a  ganaric  norm 
translation  problem  and  develops  a  set  of  tools 
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required  to r  its  solution.  Section  3  Applies 
results  of  the  previous  section  to  the  familiar  MLV 
algorithm  of  tha  Modal  Reference  typa.  Othar  algo¬ 
rithms  can  be  traatad  similarly.  Ssction  4  discus - 
sas  tha  stability  of  tha  variable  width  dead-zone 
adaptiva  systam  and,  finally,  Ssction  5  contains 
tha  concluding  remarks. 


tha  actual  transfar  function  of  a  plant  with  feed- 
hack,  dasignad  to  follow  a  referanct  modal  which  is 
prascrihad  by  the  transfer  function  M  ( s ) .  in  tha 
absence  of  unmodeled  dynamics,  its),  perfect  match¬ 
ing  is  possible. 

Whan  £(s)j* 0,  tha  stability  of  £(£,5)  can  be 
ensured  by  (requiring)  enforcing  the  condition 


2.  mathematical  PREiiMJxmies 


1+A<9,  ju)  [l+£(ju))*0 


(3) 


In  this  section  we  develop  the  necessary  tools 
for  the  definition  of  the  variable  width  dead  zone. 
As  was  already  pointed  out  in  the  introduction ,  the 
objective  is  to  find  satisfactory  bounds  for  the 
contribution  of  the  unmodeled  dynamics  to  tha  out¬ 
put  of  tha  adaptively  controlled  process,  so  that 
an  "accurate"  error  dead-zone  can  be  defined. 

The  process,  complete  with  uaaodalad  dynamics 
is  assumed  to  be  of  the  Ooyle-Stein  type,  with  the 
high  frequency  dynamics  entering  aultiplicatively; 
i.e. 

g(s)  •  g  (s) (l*£(s) ) 

P 

where  g (s)  is  the  actual  plant  transfer  function, 
g  (s)its  modeled  part  and  £(s)  tha  unmodaled  dynmn- 
ics.  Typically,  a  bound  on  tha  magnitude  of  £(s) 
is  assumed  to  be  negligible  for  frequencies  below 
crossover,  becoming  only  appreciable  for  higher 
frequencies;  no  phase  information  can  be  assumed. 
The  problem  is  then  to  find  a  bound  for  the  output 
of  the  adaptively  controlled  process  due  to  £(s). 
This  is  achieved  in  two  stages,  as  the  following 
subsections  indicate.  We  further  remark  here  that, 
due  to  feedback  in  the  adaptive  loop,  the  unmodaled 
dynamics  indirectly  influence  all  tha  state  var¬ 
iables  of  the  nominal  adaptive  loop,  with  the  mag¬ 
nitude  of  their  contribution  depending  on  the  nanare 
of  the  inputs  to  the  plant. 


2.1  Transfer  Function  Magnitude  Bounding 


In  this  subsection  we  will  derive  a  bou'  c  $ (u) 
on  the  magnitude  of  tha  frequency  response  of  a 
special  class  of  transfar  functions,  that  typically 
arise  in  MRAC  systems.  Consider  an  LTX  transfer 
function  G(9,s)  of  the  form 


G(6,s) 


MU)  {• 


i<») _ 


} 


(24) 


where 


M (s)  is  a  completely  known  stable  LT2  transit 
function 

i+Ate.sj  -  J  ts^ten  Ub) 

i-l 


<  9  <  3 
- u 


This  condition  is  satisfied  if,  for  all  9  in  tha 
space  of  admissible  parameters  tha  following  is 
true: 


|  |4(Ju)  |  |<  iQ(u)  <  l+||A"1(e.ju)||  (4) 

Assuming  this  condition  is  true,  we  proceed 
to  derive  an  upper  bound  on  |  JcIjJT  ju)  1  |  ,  para¬ 
metrized  by  |  |m ( ju)  |  j .  From  ecn7  (2), 

G(9, jw)*M( ju»>  (S) 

where  D(£,  ju)51**-A(£,  ju) 


Next,  representing  £(ju),  D(£,ju)  and  At£, jw) 
in  polar  form. 


|G(9,ju)  |  |<|  | M { ju)  ] 


1  j  £(  ju)  l|ej  * 


,0n  (£) 

l0(9.Ju)| |*J  0 

j  («,  (9KeJ 
||£(ju)A(9,  ju)  |  |e  ^ 


t6) 


An  upper  bound  fer  j  { G(£,  j.)  j  j  can  be  found  by 
maximizing  and  minimizing  respectively  the  values 
of  the  numerator  and  denominator  terms  in  (6) .  2ie 
denominator  achieves  its  smallest  value  if  the 
vectors  £(ju)A(9, ju)  and  0(9, ju)  are  oppositely 
aligned  and,  in  addition ,  —|  |  £  (ju)  1 1  achieves  its 
maximum  allowable  value  for  the  f- interval  of  in¬ 
terest.  The  above  two  conditions  axe  satisfied  if 

the  phase  angles  are  such  that  c.-C  (6)«g^(9)^^ 

*.  A  -  D  — 

and  i | £( ju) | |-£  (u) .  Mote  that  condition  (4)  for 

the  stability  of  G(£,  s)  ensures  that 
|  |£^(ju)A(9ou)  |  ,:<i  .0(9,  ju)  |  \  and,  indeed,  the 

choice  1 1  £(ju)  |  \  (u)  guarantees  minimization  of 

the  denominator  of  G (r_,  jw)  with  stability  maintained. 

Unlike  the  denominator,  the  numerator  aegnitxSe 
is  independent  of  and  is  directly  maximized  by 

choosing  | |£(ju) j j»£  U)  and  appropriately  to 

satisfy  the  denominator  phase  angle  condition. 
Consequently ,  from  eqn.  (6)  we  can  now  write  an  up¬ 
per  bound  for  G (£.  ju)  as  follows  : 


£  unknown  constant  parameter  vector  with 
~  specified  bounds  (2c) 

S>  (0)<  0  i-l . n  (2d) 

— x  **  * 

l  j£{Jw)  1 1<,  £Q(u)  for  known  i^(u)>  0  Ue) 


Note  that  when  i(s)*0,  G{£,  si«M(s),  and  therefore 
stable.  In  the  context  of  MPAS,  G(0,s)  represents 


l  <u) 
o 


i  I D(9^»  Ju)  1 t  -i  (u)  j  lA(6,jw)  ,  | 


(7) 


Finally,  we  aey  search  the  space  of  allowable  values 
of  9,  to  determine  the  desired  transfer  function 
bound 

$<u)  2  max  1 jGlgTju) i !  («) 


2 


and  eqn*.  (25)  and  (26)/  it  is  straightforward  to 
conclude  that 

£<t)«  L1  2  (27) 

From  this  point  on#  tha  proof  visas  standard 
arguments,  for  tha  boundadnass  of  u,v,£,  as  thay 
first  appaarad  in  [7]  and  outliaad  in  (4).  Ms  only 
r— »r)c  hara  that/  in  our  ca sa#  | v(t)  | <|s (t)  |  and/ 
furthermore,  v(t)-0(sup| |w(t) | | 1 ,  as  follows  from 
t>T 

eqn.  (22)  and  fig.  4.  Tha  raadar  is  again  rafarrad 
to  [6]  for  all  tha  da tails  of  tha  stability 
arguments. 
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Figure  5: 


6.  CONCLUSIONS 

A  variable  dead-xone  nonlinearity  was  introdaoad 
in  a  standard  modal  reference  adaptive  control 
algorithm  to  maintain  its  stability  in  tha  presmice 
of  unmodalad  dynamics.  Tha  variable  width  dead- 
zone  is  datarminad  on-line  on  the  basis  of  prior 
information  about  plant  parameter  bounds  and  un¬ 
modalad  dynamics  as  wall  as  about  information  ob¬ 
tained  during  adaptation.  Besides  maintaining 
stability,  the  algorithm  is  able  to  overcome  the 
conaervetiveness  of  fixed  dead -son#  on  exponential 
forgetting  factor  adaptation  mechanisms  141,  (81,  as 
simulation  results  show.  Due  to  space  limitations, 
those  are  deferred  until  tha  conference  presenta¬ 
tion  of  tha  paper. 
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In  the  following  subsection  we  will  use  $ (u)  in 
order  to  bound  the  output  of  G  (8 ,  juj)  ,  in  an  Ab¬ 
solute  value  sense,  given  an  input  x(t). 

2.2  Absolute  Value  Output  Bounding 

Consider  the  system  shorn  in  figure  1  with 
input  x  (t) ,  output  y  (t)  and  G(0_,  s)  as  defined 
before 


7  J"  ***  ii  5(£.juij,i  ;  ix(ju  I :  (ui 

But  from  the  previous  subsection  b (w)  was  detar- 
mined  such  that 

$<w)>  J |C<e, jw) | i  (15) 

Hence ,  ( 14 )  becomes , 


x(t) 


G(fts)~g  <0ft) 


y(t) 


Figure  1: 

By  taking  the  inverse  Laplace  transfer  of  G(£,s), 
we  obtain 


L'l{G(9.*)}  -  g(e.t)  (9) 

Furthermore,  since  g(6,t)  represents  a  linear-time- 
invariant  system#  by  definition  of  GCB^s)#  the 
output  y(t)  is  given  by 


y(t) 


gt^t)  • 


x(t)  • 


dx  g(8.t-x)x(x) 


<10) 


Substituting  now  in  (10)  for  g(0, t-T)  the 
expression  for  the  inverse  Laplace  transform  of 
G(£#s)  and  recalling,  further,  that  G(6.,s)  is  a 
stable  transfer  function,  we  can  write 

y(t)  ■  J  dx  — j  jj  ds  G(8,s)e#t"  T*  x(x)  • 

m  ^  t 

"hi  J  ^  a  / dT 

—  — '  Ul> 

Next,  with  u(t)  representing  the  unit  step 
function  we  have 


F  {x(t)y(-tj-  J  dx  x(T)u(-x)e*^tart»^  dx  xd)!*3^ 

(12) 

Define  F\x(t)u  (-T) }  with  the  symbol  X(ju»)  and  subs¬ 
titute  in  (11) .  Then 


y<t)  •  |r  y  dw  G<e,Jtt)3c(Jw)eiwt  U3> 


|y(t)|<  i(w)||x(ju>)||  =  y(t)  (16) 

-a» 

which  admittedly  represents  a  looser  bound  on  the 
absolute  value  of  y(t)  than  egn.  (14).  However, 
the  bound  yXt)  can  be  calculated  more  readily.  The 
inequality  (16)  has  the  following  interpr station. 
Given  a  bound  on  the  frequency  response  of  a  sys¬ 
tem,  a  bound  on  the  magnitude  of  its  output  due  to 
an  input  x(t)  can  be  calculated  at  any  instant  of 
time  by  using  the  time  history  of  the  system  input 
up  to  and  including  that  instant  of  time.  We  note 
here  again  the  time  dependence  of  T(jw)  according 
to  eqp.  (12) . 

3*  MRAC  WITH  RELAXED  TRACKING  E3P0R  Cr-TSRION 

In  this  section  we  employ  the  results  of  section 
2  to  derive  a  variable  dead -rone  width  for  the 
parameter  update  of  the  N-L-v  algorithm,  which 
ovorcommt  the  conservativeness  of  the  fixed  width 
Peterson-Narendra  scheme  [4] .  Before  we  proceed, 
we  briefly  review  the  concept  of  parameter  update 
using  a  dead-rone  nonlinearity . 

In  (4)  the  authors  have  shown  stability  of  the 
system  depicted  in  figure  2  below. 


& 


c 


v 


Figure  2 i 

In  fig.  2  the  standard  notation  is  used,  with  k 
representing  parameter  errors,  ^  filtered  auxil¬ 
iary)  state  variables  v(t)  output  deterministic 
disturbances  of  bounded  magnitude,  differentiable 
and  uniformly  continuous,  f  the  output  error  with 
disturbances!  n  represents  tbs  (part  of  the)  error 
actually  used  in  the  parameter  adaptive  law  and  is 
obtained  by  6  passed  through  the  dead- cone  of  width 
£. 

It  is  not  our  purpose  here  to  present  the 
details  of  how  the  above  error  system  as  shown  in 
fig.  2 ,  is  arrived  at.  The  reader  is  instead 
referred  to  (4)  for  those  as  well* as  the  stability 
proof  of  that  modified  algorithm.  We  simply  present 
here  the  parameter  adaptive  lavs,  with  tne  dead- 
zone  nonlinearity  in  their  simplest  form,  for  the 
sake  of  completing  the  proo*»m  iescriptior.  v.-uer. 
forms  the  basis  for  the  developments  ir.  the  present 


By  the  Cauchy Schwartz  inequality  it  follows  from 
(13)  that 


paper.  The  ptfimur  adjustment  is  as  described 
by  eqns.  (17)  below. 


Jett)  -  E(t)  - 


with  Hit) 


axj 

l-K»CT«t)Clt) 

«(t)  )e<«)|>e 

o  |e(t)|<E 


(17a) 


(17b) 


where  £  represents  the  width  of  the  deed-zero. 

Although  in  [41  the  discussion  is  not  perticulerly 
enlightening  as  to  how  exactly  the  magnitude  £  is 
decided  upon,  it  is  the  present  Authors' 
conclusion  that  most  likely,  in  [4] 

E  >  I |v(t) I |„  «  B«x|v(t) I  (18) 

t 

end,  therefore,  is  very  conservative  as  the  seat 
authors  have  pointed  out  in  [51 .  we  next  proceed 
to  analyse  the  original  NLV  Model  Reference  algo¬ 
rithm,  as  represented  in  Fig.  3,  with  the  plant 
dynamics  now  replaced  by  the  actual  plant 
P(6,s)  tl+£(s)]. 


for  P(9,s) .  That  is,  the  designer  knows 

(i)  an  upper  bound  on  the  relative  degree  n*  of 
P(9,s) 

(ii)  an  upper  bound  on  the  degree  n  of  F  (€,,*/ 

(iii)  that  P^,*)  is  minimum  phase 

(iv)  the  sign  of  the  high  frequency  gain  of  ?(£,*) 

the  £(s)  part  of  the  plant  represents  the  (multi¬ 
plicative)  uncertainty  associated  with  the  nominal 
plant  P(0,s).  This  uncertainty  is  due  to  high  fre¬ 
quency  dynamics,  which  are  assumed  of  unspecified 
structure  but  satisfy  a  magnitude  constraint 

||t(jw)||<  l  (ui) ,  as  already  mantioned.  Further, 

*  o 

we  note  that  the  "actual"  plant  representation 
P(9,s)  (l+i(s)J  will  not  in  general  satisfy  any  of 
the  four  standard  assuttptions  listed  above.  This 
fact  becomes  pivotal  in  the  inability  of  the  adap¬ 
tive  controller  to  achieve  transfer  function  match¬ 
ing  between  the  compensated  plant  and  the  reference 
model.  As  a  result,  it  becomes  impossible  for 
general  inputs  to  drive  the  tracking  (output)  er¬ 
ror  to  zero.  However,  one  may  expect  the  tracking 
error  to  be  small,  if  the  plant  is  excited  by  sig¬ 
nals  with  dominant  low  frequency  content  over  the 
range  where  ?(£,*}  is  a  good  approximation  to  the 
actual  plant  transfer  function. 

In  what  follows  we  will  next  show  that  the 
error  system  underlying  the  structure  in  figure  3 
differs  from  that  of  the  same  structure,  as  shown 
in  fig.  2,  where  £(s)»0,  only  by  an  additive  per¬ 
turbation  ears  in  the  output.  This  term  can  be 
bounded  using  the  results  of  Section  2  and  a  var¬ 
iable  width  dead-zene  can  be  defined  for  the  adap¬ 
tation  mechanism.  Stability  of  the  scheme  is 
subsequently  discussed  in  Section  5. 

We  start  by  considering  first  the  case  where 
£(s)*0  in  Figure  3..  In  this  case  the  standard 
MRAC  assumptions  about  the  plant  are  true.  There 
exists  a  vector  k»(6)  of  fixed  gains  which,  when 
applied  to  the  system,  results  in  matching  of  the 
compensated  plant  transfer  function  with  that  of 
the  model .  We  adopt  the  shorthand  notation  (6J 

and  C.  (9J  to  indicate  the  LTI  transfer  functions 
(k£(9.) )  and  (k*  (£)  ,k*  (9) )  respectively.  Now 

assuming  i(s)^0  but  maintaining  the  same  defini¬ 
tions  of  Cx(k*(e)  »  Cx(9)  and  C2  (k*  (e)  ,k*  (9) )  ,  • 

Cj ( £)  based  on  the  reduced  model,  we  may  derive 

an  expression  for  the  error  system  as  follows, 
where  for  convenience,  the  argument's'  has  been 
suppressed  throughout. 


Figure  3: 


The  component  £,(9,*)  incorporates  the  designer's 
knowledge  of  the  dominant,  low  frequency  response 
of  the  plant,  including  a  vector  9  of  uncertain 
parameters,  known  only  within  crecooputamle  bounds. 
The  standard  MRAC  assumptions  aoout  the  plant  hold 


V  C, (0)P(9) [1*£] 
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R  I+cTlafcI  (6)P<6)  [1+i.J 

V! >f(8>  cx<e)p<s) 

i-*c,  (e)c,(8)p(9) *  I^r<iTc^irp(ir 
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1«-C.  (9)C2(9)P(9)  il+i]  (19) 

3y  definition  of  C. (z)  and  C, (9)  we  have 

1  “  * 

C.  (9)P<6) 

(.♦C,  <i)C,(S)P(S)  fl  :2C1 

i  ^  *  *  “ 


(24) 


Using  this  fact  and  introducing  the  notation 
A(0)  2  C^C (6.)? (£)  «•  can  writ*  eqn.  (19)  in  a 

more  compact  fora  aa 


E(t) 


5  J  t  (U)  | ! 


U(ju)  | 


2n-i 

I 

i»l 


!  vT  t  juj> 


%-**“•  langfeu —  (2l) 

Next#  defining  £(t)  *  k(t)-k*#  referring  to  fig.  3 
and  interchanging  time  rtnmain  and  tranaforaad 
quantities,  we  dariva  an  expression  for  €  aa  givan 
in  eqn.  (22) . 

yB  yn  .  y_ 

€-«♦  +  *ETw--£  •  R 

R  R  —  R 

T  IT  U2) 

a  Z-l*  *  '  l+A  (6)  (1+iJ  '  (-  -*r) 

For  tha  caaa  where  Ms>-0#  this  rasult  raducas  to 
the  standard  augnantad  HJRAC  arror  system  of 
Narendra,  Lin  and  Valavani  with  L*1  *  N.  Tha  new 
arror  system  is  shown  in  Figura  (4)  with  a  variable 
dead  zona  non* linearity  added  to  tha  output  signal 
path. 


Figure  4: 


Ha  observe  that  tha  system  is  of  tha  form  shown  in 
fig.  2.  In  order  to  specify  a  stable  adaptive 
law,  wa  need  to  find  a  bounding  signal  E(t)  >j  v(t)|Yt. 
Ha  may  redraw  tha  arror  system  using  the 

fact  that  £*k#-k  and  tha  input  to  tha  plant  in 
*  “  *“■ 

Figura  5  is  u*k  w  +r.  The  resulting  represen¬ 
tation  is  shown  IrTpigure  5  below, 
rrom  precomputed  bounds  on  8,  bounds  on  k*  can  be 
precomputed  also.  He  now  make  the  definition 


k*  ■  | max  k* 
8  ~ 


(23) 


whara  tha  maximisation  over  9^  is  carriad  out  in¬ 
dividually  over  every  component  of  k*.  Using  (23) 
in  conduction  with  tha  results  of  Section  2  it 
readily  follows  that  an  upper  bound  £(t)  for  v(t> 
can  be  computed.  More  specifically,  wa  can  write 


whara  eqn.  (12)  has  bean  used  for  calculation  of 
tha  transforms  ITtjw)  and  STtjw)  for  tha  input  u(t) 

and  signals  u^tt)  respectively.  An  adaptation  law 

of  tha  fora  described  in  egn.  (17)  can  than  be 
employed  with  tha  width  of  tha  dead  zona  defined 
by  eqn.  (24) .  The  resulting  schema  is  stable  and 
an  outline  of  its  stability  proof  is  given  in  the 
following  section. 


5.  STABILITY 


The  stability  proof  of  the  proposed  algorithm 
with  variable  dead-zone  follows  along  very  similar 
lines  for  the  most  part  with  that  in  [4]  .  However, 
in  the  present  case  it  is  additionally  conditioned 
on  the  reference  model  definition  and  the  admissible 
parameter  set,  as  egn.  (4)  of  section  2.1  implies. 
More  specifically,  the  space  of  admissible  para¬ 
meters  is  implicitly  defined  through  the  reference 
model  by  eqn.  (4),  in  conjunction  with  condition 
(2c)  and  is  such  that  the  desired  (class  of)  ref¬ 
erence  model (s)  remains  stable  in  the  presence  of 
the  unmodeled  dynamics  Ms)  of  the  plant.  This  is 
a  standard  and  reasonable  assumption  made  in  the 
design  of  all  fixed  parameter  controllers  as  well. 
Due  to  space  considerations  we  will  not  elaborate 
on  this  further  but  will  instead  refer  the  reader 
to  [6]  for  store  details. 

Consequently,  given  eqn.  (4) ,  which  is  fun¬ 
damental  tvan  for  a  non- adapts vs  design,  the  effect 
of  unaodeled  dynamics  can  be  represented  as  an  out¬ 
put  perturbation  v(t)  as  suggested  in  ecr..  (22)  and 
depicted  in  fig.  4.  v(t)  is  the  output  of  a  stable 
linear  system  which  is  bounded  for  bounded  inputs. 
He  next  proceed  to  outline  the  steps  for  proving 
boundedness  of  *,u,;,v  and  the  output  arror  €. 

The  boundedness  of  £  follows  directly  from 
the  standard  Lyapunov  function  dafinition 

(V(k)  -  j  (£*£)  and  the  adaptation  law  ,17a)  in 

conjunction  with  eqn.  (22)  where  E(t)  is  defined. 
Also,  from  the  definition  of  the  dead- zone,  the 

term  £  k  can  be  bounded  above  and  below  by  bounds 
of  the  form 


ce (t)  J  |ST£*w(t)  f. I£(t) !  •  ;5*i*v(e> ! 

(25) 

where  f,  (•)  and  f.  (• )  are  continuous  functions 

j.  2 


of  E (t) • 

Next,  by  the  definition  of  the  Lyapunov  func¬ 
tion,  its  time  derivative,  in  conjunction  with  the 
adaptation  law  given  by  eqn.  (12a) ,  can  be  written 

as  « 


V.  E*S 


(26) 


From  the  fact  that 


/ 


V(T)dT  <• 


